Intermediate temperature superfluidity in an atomic Fermi gas with population imbalance 
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We derive the underlying finite temperature theory which describes Fermi gas superfluidity with population 
imbalance in a homogeneous system. We compute the pair formation temperature and superfluid transition 
temperature T c and superfluid density in a manner consistent with the standard ground state equations, and 
thereby present a complete phase diagram. Finite temperature stabilizes superfluidity, as manifested by two 
solutions for T c , or by low T instabilities. At unitarity the polarized state is an "intermediate temperature 
superfluid". 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

Excitement in the field of ultracold Fermi gases has to do 
with their remarkable tunability. As a magnetic field is var- 
ied from weak to strong this system undergoes a transition 
from a Bose Einstein condensate (BEC) to a BCS-based su- 
perfluid. Recently another tunability has emerged; one 
can vary the population imbalance between the two "spin" 
species. This capability has led to speculations about new 
phases of superfluidity, quantum critical points, and has reper- 
cussions as well in other sub-fields of physics yl Hi. A nice 
body of theoretical work on this subject 0, HI has fo- 
cused on zero temperature (T) studies of the simplest mean 
field wavefunction 1 9], with population imbalance. Additional 
important work presents 1 10] a T = two-channel, mean- 
field approach for very narrow Feshbach resonances, as well 
as a study of finite T effects 1 1 1], albeit without a determina- 
tion of superfluid order. 

Superfluidity is, generally, a finite T phenomenon and it 
is the purpose of the present paper to explore finite temper- 
ature effects HI Q based on the BCS-Leggett ground 
state with population imbalance. We determine the behavior 
of the pair formation temperature T*, the superfluid transition 
temperature, T c , and superfluid density n s (T), and, thereby 
arrive at a phase diagram which addresses general T. Im- 
portantly, we find that in the fermionic regime, superfluidity 
exists at finite T (although not at T = 0), leading to the new 
concept of an "intermediate temperature superfluid". Because 
temperature acts in this rather unexpected fashion, we reduce 
the complexity and confine our attention to the homogeneous 
system. 

The approach which we outline below, importantly, in- 
cludes what we call "pseudogap effects". For T / 0, the 
excitation gap A is different from the order parameter, due to 
the contribution to A from noncondensed pairs 1121 llal . In 
this way, the solution for T c is necessarily different from that 
obtained in the literature. Generally, A 2 contains two additive 
contributions fl3lfl4ll from the condensate ( A 2 C ) and noncon- 
densed pairs (Ap_), and they are proportional to the total, con- 
densed, and noncondensed pair densities, respectively. This 
decomposition is analogous to the particle number constraint 
in ideal BEC. We emphasize that the central equations derived 
below are not compatible 1 15] with the T ^ formalism of 
Ref. II 611. In addition, the "naive mean field theory" with the 
unphysical assumption that A(T) = A SC (T) is not a correct 
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rendition of T ^ effects associated with the BCS-Leggett 
ground state. 

We define the noncondensed pair propagator, as t(Q) = 
U/[l + V~x{Q)], where, as in Ref. 11411 . the pair suscep- 
tibility, given by X (Q) = ^xi^iQ - ^G^K) + 
Gqi(Q — K)G^(K)], can be derived from equations for the 
Green's functions , consistent with the BCS-Leggett ground 
state equations. Here G a (K) and Go. a (K) = iuj n — £k.<r 
are the full and bare Green's functions (with a =f, |, £k,<r = 
£k — Mct)- We adopt a one-channel approach since the 6 Li 
resonances studied thus far are broad and consider a Fermi 
gas of two spin species with kinetic energy ek = h 2 k 2 /2m 
and chemical potential ^ and p^, subject to an attractive con- 
tact potential (U < 0) between the different spin states. We 
take h = 1, k,B = 1, and K = (iw n ,k), Q = (if2 n ,q), 
J2k = ^EnSk' etc ' where u> n (Q, n ) is the standard odd 
(even) Matsubara frequency. 

In the superfluid state, the "gap equation" is given by 
U^ 1 + x(0) = 0, which is equivalent to fj, pa i r = 0, the 
BEC condition of the pairs. Below T c the self-energy can 
be well approximated |12] by the BCS form, T, a (K) = 
-A 2 G ,a(-K), where a = -a. Therefore, G hl (K) = 



■J{iuj n ±h - E k ) 



iy(ico n ± h + Ek), where E k — 

y/& + A2 ' & = £k-M,M= Or+MlJA h = Or-MlJA 
and itk'^k = (li £k/-Ek)/2. Since the polarization p > 0, 
we always have h > 0. 

The "gap equation" can then be rewritten in terms of the 
two-body s-wave scattering length a, leading to 
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where /(x) = [/(x+h) + /(a;-ft)]/2,andm/47ra = 1/U + 
^ k (2ek) _1 - Here f(x) is the Fermi distribution function. 
We define pn = Sn = n-\ — > 0, where 71 = + raj 
is the total atomic density, and p = Sn/n. Similarly, using 
n a = J^k ^cr(^), one can write 
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Note that, except for the number difference [Eq. (!2bH . all 
equations including those below can be obtained from their 
unpolarized counterparts by replacing f(x) and its derivative 
f'(x) with f{x) and f'(x), respectively. 

While Eqs. Q-(|2) have been written down in the literature 
IHIUI]- the present derivation can be used to go further and to 
determine the dispersion relation and the number density for 
noncondensed pairs. We find 



t(Q) 



which vanishes at T = 0, where A 2 
phase, t-\Q) = x(Q) ~ x(0) « 



= A 2 

z(n - 
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In the superfluid 



der in il, and after analytical continuation (if2„ 



to first or- 
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Here X (Q) - £ k [ B k +« q _ k -in„ "k ~ s k -£,_ fc +in„ "k 
It follows that the inverse residue Z = [n-2 /(£ k )]/2A 2 . 
Thus A 2 g = Z -1 6(f2 q ), where 6(x) is the Bose dis- 
tribution function. To lowest order in q, the pair dispersion 
= Q 2 /2M*, where the effective pair mass M* can be 
computed from a low q expansion of £! q . This q 2 dispersion 
is associated 11211 with BCS-type ground states, which have 
been the basis for essentially all population imbalance work. 

Importantly, Eqs. Q-Q can be used to determine T c as the 
extremal temperature(s) in the normal state at which noncon- 
densed pairs exhaust the total weight of A 2 so that A 2 g = A 2 . 
Solving for the "transition temperature" in the absence of 
pseudogap effects leads to the quantity T C MF . More precisely, 
T^ F is defined to be the temperature at which A(T) vanishes 
within Eqs. Q and (0- This provides a reasonable estimate 
for the pairing onset temperature T*, when a stable super- 
fluid phase exists. It should be noted that T* represents a 
smooth crossover rather than an abrupt phase transition, and 
that Eq. Q must be altered |17] above T c to include finite 
Impair ■ We will see that understanding the behavior of T^ F 
is a necessary first step en route to understanding the behavior 
of T c itself. 

The superfluid density n s (T) is also required to vanish at 
the same value(s) for T c , as deduced above. Our calculation 
of n s closely follows previous work 1141 1 1 811 for the case of 
the unpolarized superfluid. There is an important cancellation 
between the current vertex and self-energy contributions in- 
volving A 2 g so that, as expected, n s (T) varies with the order 
parameter A 2 C . It is given by 
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which at T = agrees with Ref. j5[. 

The stability requirements for the superfluid phase have 
been discussed in the literature Q. In general, one requires 
that the superfluid density be positive and that the 2x2 "num- 
ber susceptibility" matrix for dn a /dp, a i have only positive 
eigenvalues when the gap equation is satisfied. The A de- 
pendence of n a introduces into the matrix the overall factor 

( Is' ) ■ Thus, the second stability requirement is equiva- 
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Figure 1 : (Color online) Mean-field behavior of T c as a function 
of 1/Uf a for different p. Shown in the inset is the pairing gap A (T) 
at different l/kpa for p = 0.3. Here Ef = fesTp = h 2 k F /2m is 
the noninteracting Fermi energy for p — 0. 



lent to the condition that 
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Here fl is the thermodynamical potential. A third stability 
requirement, specific to the present calculations, is that the 
pair mass M* > 0. 

In Fig.Q we present a plot of T^ IF as a function of 
for a range of p. In the inset we plot A(T) at different 1/kpa 
forp = 0.3. Forp < 0.9 and sufficiently low T C MF , the curves 
for T C MF develop an unexpected structure, as one sweeps to- 
ward the BCS regime. Once 1/kpa is less than a critical 
value, (l/fcpa) c where T C MF vanishes, there are two T^ F 
lines. The lower branch starts from (l/fcj?a) c and increases as 
1/kpa decreases. This structure implies that A is nonmono- 
tonic 1 19] in T, as indicated by the bottom curve in the inset 
of Fig. ^ The two zeroes of A represent the two values of 
T C MF . In contrast to the more conventional behavior (shown 
in the top curve for stronger pairing interaction), A increases 
with T at low temperature when 1/kpa is sufficiently small. 
This indicates that temperature enables pairing. This was also 
inferred in Ref. In general superfluids, one would argue 
that these two effects compete. 

Insight into this important phenomenon in the fermionic 
regime (// > 0), is provided by studying the momentum dis- 
tribution n a (k) at T = and finite T. At T = 0, pairing is 
present only for ek below e± = Max(0, /x — \Jh 2 — A 2 ) and 
above €2 = p, + Vh 2 — A 2 . This polarized T = state re- 
quires that pairs persist to relatively high energies ek > £2, as 
a result of the Pauli principle which pushes these states out of 
the "normal" regime occupied by the majority species. This 
kinetic energy cost competes with the gain from condensa- 
tion energy and for sufficiently weak attraction this "breached 
pair" structure 1 20] becomes unstable at T = 0. By contrast, 
at finite T the regime originally occupied exclusively by the 
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Figure 2: (Color online) T c as a function of 1/kpa for different p. 
The T c curve splits into two disconnected curves for 0.14 < p < 
0.185. The T c solution inside the shaded area is unstable. 



majority species between t\ and 62 is no longer completely 
filled and pairs can "spill over" from both lower and higher 
energy states into this regime. This not only helps lower the 
kinetic energy but allows the "normal" regime to participate 
in pairing and thus lowers the potential energy. In this way 
temperature can enhance pairing. It should be noted that the 
majority species between ei and £2 contains a pairing self- 
energy and is different from a free Fermi gas. 

Figure[2]represents solutions for T c of our central equation 
set [Eqs. Q-(|3j] as a function of 1/kpa for the entire range 
of p. If the solution for T c falls into the shaded region, there is 
no stable superfluid (since d 2 il/0A 2 < 0, through Eq. (|5j). 
For low polarizations p < 0.185, the behavior of T c is similar 
to thatofT c Mi " when one approaches the BCS regime. There 
may be one or two T c 's which, when stable, will be associated 
with intermediate temperature superfluidity. Whenp > 0.185, 
however, no solution can be found for the regime 1/kpa < 
0.18, because M* < there. We stress that the origin of 
the intermediate temperature superfluid we find here lies in a 
very early stage of the calculations; it can already be seen as 
a consequence of the constraints imposed on the pairing gap 
in the low T regime when there is a delicate energetic balance 
between normal and paired states [see, Eqs. and and, 
Fig.ID. 

In Fig.|5]we summarize our observations in the form of a 
general temperature phase diagram. In region I, the system 
is normal and superfluidity is absent. However, this normal 
phase need not be a Fermi gas. Close to the boundary, as 
shown in the inset to Fig. [2 (bottom curve) finite T pairing 
may occur with or without phase coherence. In region III, sta- 
ble superfluidity is present for all T < T c . Finally in region 
II, within the shaded region (IIC and IID), we find a stable po- 
larized superfluid phase for intermediate temperatures, not in- 
cluding T = 0, which we refer to as intermediate temperature 
superfluidity. In IIA and IIB no stable polarized superfluid is 
found. The nearly vertical blue line shown in the figure rep- 
resents the line T c = which appears around 1/kpa 0.18 
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Figure 3: (Color online) Phase diagram on the p- 1/kpa plane with 
nearly vertical line corresponding to T c = 0. Yellow region (shaded) 
corresponds to intermediate temperature superfluidity. Region I the 
system is normal; IIA, M* < so no solution for T c exists; IIB: 
Solution for T c exists but superfluid state is unstable; IIC: Superfluid 
state exists at intermediate T between upper and lower T c 's but not 
T = 0; IID: Superfluid state exists at finite T but becomes unstable 
at a low temperature T unsta bi e (shown in the inset for 1/kpa = 0.5) 
where n s is finite; III: Superfluid state exists for all T < T c > 
0. The chemical potential fi changes sign within regions IIB and 
IID (close to 1/kpa ~ 0.6 for all p), while n s (0) vanishes along a 
nearly vertical line between (p, 1/kpa) — (1, 0.3) to (0, 0.6). The 
p = boundary is not continuously connected to the rest of the phase 
diagram. 



(See Fig. |3 and is roughly independent of polarization. 

Finite momentum condensates may well occur in 

any of the regimes in II, particularly IIA and IIB for which 
our equations do not yield stable zero momentum condensa- 
tion. Future work will explore the nature of the stable phases 
in these regimes. The boundaries of the region denoted II can 
be compared with other T = phase diagrams in the liter- 
ature la. [Toll . In contrast to Ref. @L we find that the most 
stringent criterion for stability at T = is the positivity of 
the second order partial derivative of d 2 il/dA 2 [as given in 
Eq. @], This defines the boundary between II and III. This is 
substantially different from the line associated with n s (0) = 
(used in Ref. |5]) which is described by a nearly vertical line 
from (p, 1/kp a) = (1, 0.3) to (0, 0.6). Similarly the locus 
of points in the two-dimensional parameter space (p, 1/kpa) 
where T^ F vanishes defines the boundary between I and II, 
as is consistent with its counterpart in Ref. |g]. 

In Region IID, we define T unsta bie (which is below the sin- 
gle T c ) as the temperature where the system becomes unsta- 
ble, via Eq. Q. At a given 1/kpa, T unsta bie decreases with 
decreasing p, and approaches as p — > 0. This is shown in the 
inset to Fig. 3 for 1/kpa = 0.5. In region IIC, the lower T c 
approaches as p approaches 0. Thus, there is an important 
distinction between the p = andp — * + limits, especially at 
T = 0. For p small but finite, calculations readily encounter 
instabilities at strictly T = and here superfluidity is very 
fragile to the introduction of small imbalance. By contrast, at 
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Figure 4: (Color online) Normalized superfluid density n a /n as a 
function of T/T F atp = 0.1 for various l/k F a from BCS to BEC, 
corresponding to regions IIC {l/kwa = —0.5 and 0), IID (0.2, 0.5, 
1.0) and III (1.5), respectively. The inset plots n s o = n s (T — 0) 
versus 1 — p for different 1/kpa, indicating n s o — > 2n^ in the BEC 
limit. The dotted (segments of the) curves represent unstable solu- 
tions within region IIB (IID). 

finite T this fragility is not as pronounced. We conclude that 
only T = is a problematic temperature for weakly polarized 
superfluidity. 

Figure 0]presents n s (T) for p = 0.1. Except for the case 
1/kpa = 1.5, Fig. |4] shows the typical behavior in region 
II (of Fig. [3}, corresponding to intermediate temperature su- 
perfluidity. The observations here (and associated nonmono- 
tonicities) for n s (T) are similar in many ways to what is seen 
for A(T) in the inset to Fig. [2 In region IIC, n s goes to zero at 



the upper and lower T c , whereas in IID, n s (T) abruptly stops 
at T unsta bie- The dotted lines indicate that they are in the un- 
stable regime. Throughout region III, n s is found to be mono- 
tonically decreasing with increasing T, as in conventional su- 
perfluids. Finally in the inset of Fig.|4]we plot n s o = n s (0) as 
a function of 1 — p = 2n j /n. Only in the deep BEC regime 
is the dependence linear. This plot reflects that the excess un- 
paired fermions interact with the paired states, leading to a 
reduced superfluid density at T = relative to 2nj,. 

The experimental situation regarding the stability of a uni- 
tarypolarized superfluid (UPS) is currently being unraveled 
iQyE- If one includes the trap, within the local density ap- 
proximation it appears |0, 0, [Till that the local polarization 
p(r), in effect, increases continuously from a small value at 
the trap center to 100% at the trap edge. It follows from this 
paper that at very low T the superfluid trap center will not 
support polarization, but for a range of T closer to T c , po- 
larization can penetrate the core. We estimate from Fig. [2] 
(assuming the central p w 0.05), that there exists a UPS for 
T ~ 0.05 — 0.25X/.'. Given the temperature range in experi- 
ment J]]] this appears to be not inconsistent with current data 
(Tp = 1.9/zK, T = 300 ~ 505nK = 0.16 - 0.27T F on res- 
onance). In the near-BEC regime our predictions also appear 
consistent with new data in Ref. [22]. More generally, because 
the local p(r = 0) is small the unstable region is suppressed 
to very low T as p — > 0, this may explain why superfluidity 
in atomic traps can be observed experimentally. Future the- 
ory including the trap will be required to provide quantitative 
comparison with experiment. 
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